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ABSTRACT 

Theoretical work has shown that intermediate mass (O.OIMq < Mue < O.IMq) Helium shells will unstably 
ignite on the accreting white dwarf (WD) in an AM CVn binary. For more massive (M > O.8M0) WDs, these 
helium shells can be dense enough (> 5 x 10^ g cm"^) that the convectively burning region runs away on a 
timescale comparable to the sound travel time across the shell; raising the possibility for an explosive outcome 
rather than an Eddington limited helium novae. The nature of the explosion (i.e. deflagration or detonation) 
remains ambiguous, is certainly density dependent, and likely breaks spherical symmetry. In the case of det- 
onation, this causes a laterally propagating front whose properties in these geometrically thin and low density 
shells we begin to study here. Our calculations show that the radial expansion time of < 0.1 s leads to in- 
complete helium burning, in agreement with recent work by Sim and collaborators, but that the nuclear energy 
released is still adequate to realize a self-sustaining laterally propagating detonation. These detonations are 
slower than the Chapman- Jouguet speed of 1.5 x 10'^ cm s"', but still fast enough at 0.9 x 10"^ cm s~' to go 
around the star prior to the transit through the star of the inwardly propagating weak shock. Our simulations 
resolve the subsonic region behind the reaction front in the detonation wave. The 2D nucleosynthesis is shown 
to be consistent with a truncated ID Zeldovich-von Neumann-Doring (ZND) calculation at the slower detona- 
tion speed. The ashes from the lateral detonation are typically He rich, and consist of predominantly '^'^Ti, "^^Cr, 
along with a small amount of ''^Fe, with very little ^^Ni and with significant '*"Ca in carbon-enriched layers. If 
this helium detonation results in a Type la Supernova, its spectral signatures would appear for the first few days 
after explosion. 

Subject headings: hydrodynamics — nuclear reactions, nucleosynthesis, abundances — white dwarfs 



1. INTRODUCTION 



Re cent theoretical insights ([Bildsten et al.|2007[|Shen et al. 

IWl 



20T0l[Waldman et al.|2010[ [Woosley & Kasen|2011| i and ob 
servations of ve ry rapidly evolving therm onuclear supernovae 
such a s 2002bj (Poznanski et al. 2010), SN2010X fesliwal] 
[eFaL][2010 ), SN 1939B, and SN 1885A (Perets et al...20l]T 
motivate an investigation of the outcome of unstable burn- 
ing of intermediate mass (O.OIM© < Mhb < O.2M0) Hehum 
shells on accreting white dwarfs (WDs). |Shen & Bildsten 



(|2009 ) showed that the burning initially proceeds as a spheri- 
cally symmetric slo wly evolving (minutes to ho urs) convec- 
tive layer [see also [Woosley & Kasen] ( |201 l| l]. For even 
lower shell masses (such as occur in stars on the Asymptotic 
Giant Branch), the convective burning causes radial expan- 
sion adequate to increase the evolutionary timescales to times 
long enough for confidence in a one-dimensional hydrostatic 
calculation. However, there is a critical helium shell mass 
(Mhs > 0.1(0.01)Mq for a WD mass of M = 0.6(1. 1)Mq) 
above which the convective burning timescales become so 



short that a dynamical mode of burning will occur ( Shen & 

[Bildsten 2009). 

As highlighted by many (Bildsten et al. 2007; She n et al.| 
|2010[ [Woosley & Kasen 201 1 ), the dynamically burning He 
shell differs from the well studied case of dynamic carbon 
burning in WD cores during Type la SNe, as the finite gravi- 
tational acceleration, g, in the He shell rapidly forces mixing 
of any unstable density gradients from a radially propagat- 
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ing deflagration ( [Timmes & Niemeyer|2000"| ). This prohibits 
the localization of a well defined and propagating deflagra- 
tion front. If detonation is the result, we expect that a lack of 
synchronicity across the star will lead to laterall y propagat- 
ing detonation wave(s) from the ignition site(s). Fink et al. 
p007, 2010); Sim et al. (2012) revived the exciting po ssibil- 
ity (|Liv ne & Glasner 1991 ; Woosley & Weaver 1994; Livne| 
& Arnett,,1995) that the weak shock that penetrates the un- 
erlying WD from a laterally propagating He detonation will 
strengthen enough (even for a Mne = O.OIM0 shell) at the 
off-center focus point to trigger a C/0 detonation in the WD. 
They have argued that nearly all He detonations will lead to 
a complete explosion of sub-Chandrasekhar C/O mass WDs. 
This motivated theoretical modeUng of such explosive events 
(|Sim et al.||20T0l [Kromer et al.|[20T0l [Waldman et al.|[20T0l 



Woosley & Kasen|20 1 1 Sim et al.l2()12 1 finding similarities 
to the width-luminosity relation or Type la SNe and poten- 
tially offerin g new insights on u nusual events [e.g. S N 2005E 
( |Perets etal.|201 0), PTF09dav ( |SuUivan et aLpOll) , and odi- 
ers dKasUwal et al. 2011)1 

Though un certainties remain ( [Shen et al.|2010[|Woosley & 
|Kasen|[T01 1 } as to the exact criteria that rnust be met to ini- 
tiate a detonation in these lower density (p < 10^ g cm"^) He 
shells, we work here to resolve how the lateral propagation 
will be impacted by the radial expansion that occurs during 
the burning. The radial blow-out timescale is H/cs ^ 0.01 s 
for a shell of thickness // « 10^ cm, where the sound speed, 
Cs ~ 10^ cm s~\ is set by the post-shock temperature of the 
detonation front. This timescale decreases for thinner layers, 
whereas the nuclear burning timescale in the detonation reac- 
tion zone increases. If a propagating detonation wave remains 
stable under incomplete burning, it would reduce the propa- 
gation speed to values lower than the Chapman-Jouguet (CJ) 
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speed of 1. 5 x 10^ cm s ' c alculated for a complete burn of 
He to ^^Ni ( Sim et al. 2012 1. At even lower densities (or shell 
thicknesses), the blow-out occurs so rapidly that the propaga- 
tion could be quenched. 

We begin in ^by showing a 2D simulation of a steadily 
propagating detonation in a thin helium layer for an example 
case motivated by the AM CVn accretion scenario. In or- 
der to investigate the robustness of our numerical techniques 
and better understand the incomplete burning that occurs, we 
undertake two sets of comparisons. In ^ we perform ID 
and 2D simulations, with enough resolution to resolve all 
the burning scales, investigating transient strengthening det- 
onations in uniform pure He. We find that the incompletely 
burned strengthening detonation provides a useful compari- 
son to the blowout case. In ^ we compare our simulations to 
integrations that give the structure of ID steady-state detona- 
tions (ZNEq performed with a large (200-nuclide) reaction 
network ana find qualitative agreement. Finally, in ^ we 
discuss astrophysical implications and future work. 

2. SURFACE BLOWOUT INDUCES FREEZE-OUT OF 
INTERMEDIATE PRODUCTS 

Our fiducial model for the multidimensional runs was mo- 



tivated by a tirn e dependent M calculation with MESi^tPax- 
|ton et al.ll201l!) that had a declining M as expected from a 



^generate He donor in an AM CVn binary (Bildsten et al. 

calculation had! 10 He Hashes 



|2006j|2007|. This exploratory 
at very high accretion rates o 



very iiigli accretion rates on the M = 1 .QMq WD, with 
about two-thirds of the accreted material leaving during the 
He novae. The remaining carbon-oxygen ashes from the He 
burning creates a hot (7 sa 8 x 10^ K) layer of sa O.OISM©. 
The final flash occurred when M « 6 x 1O~**M0 yr"', with a 
freshly accumulated helium layer mass of O.O28M0. At this 
M, the temperature maximum occurs above the base of the 
He shell, placing the base of the convective zone at a mass 
coordinate inward from the surface of O.O14M0. Our fidu- 
cial model is co mparable to two M = I.OMq WD models of 
IWoosley & Kase n (2011 ). The first (lOAA) had a cold C/O 
core (r < 3 X lO"^ K) and M = 7 x IO'^^Mq yr"\ giving an ac- 
cumulated mass of M„ = Q.Q52Mq and a convective layer of 
O.O262M0, thicker than our He layer due to the colder core. 
The second (lOHB) had a hot C/O core {T « (7-8) x 10^ K) 
and M = 5 X 10"** yr"', giving an accumulated mass of 
Ma = O.O22M0 and a convective layer of O.O22M0, much 
closer to our fiducial. Compared to the models shown in their 
Figure 2, our model has a peak luminosity of lO'**^ erg s"' 
when the conditions at the base of the convection zone are 
= 0.14 and Tg = 5.0. This places our model above their 
threshold for shell detonations. 

For our simulation, we used the temperature and density 
structure when the timescale for temperature change in the 
convective burning layer was at its minimum, about 5 sec- 
onds. This calculation is intended to demonstrate one case in 
which the burning scales are long enough for blowout to trun- 
cate the reactions. This radial (vertical) structure is shown in 
Figure [T] The initial helium layer is in hydrostatic balance, 
with thenot (Tg = 0.52), carbon-enriched iXi2 = 0.1 due to He 
burning) convection zone above the colder (Tg = 0.08) pure 

* Throughout we use Zeldovich-von Neumann-Doring (ZND) to refer to 
the eigenvalue method for computing the spatial structure of steady state det- 
onations. See, e.g., Khokhlov 1 19891 for a convenient form of the equations 
integrated and Fickett & DaviS|^1979^ for a broader discussion. 

^ http://mesa.sourceforge.net 
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He layer. The base of the cold He layer defines the zero of 
the vertical coordinate, and has a density of ps = 5.89. The 
cold layer extends up to sa 0.25 x 10** cm, and has a col- 
umn depth of lO'^* g cm"^ or « O.OIMq, while the convec- 
tive layer extends to sa 1.4 x 10** cm and has a similar mass. 
Above this is a low-density "fluff layer outside the star ex- 
tending to the edge of the grid. We work in plane parallel, 
with g = 7.32 X 10** cm s"^, appropriate to the IM© WD of 
radius R = 4.3 x 10^ cm. The helium is sitting on a presumed 
(for this simulation) inert core. 

Reactive hvdrodvnarnics si mulations were performed using 
FLASI^( Fryxell et a l.'2600^ and utilized the "approxl3" nu- 
clear reaction network (Timmes et al. 2000a) that accounts 
for the fas t (a,p)(p,j ) side channels relevant to alpha burn- 
ing (Timmes et al.| 2000b). We used the pi ecewise parabolic 
method (PPM; Colel la & Woodward|1984| l for compressible 
hydrodynamics. Reactions w ere disabled in regi ons of the 
grid within a detected shock (Fryxell et al. 1989 1. The hy- 
drodynamic and reaction evolution are operator split and the 
temperature is held fixed during the reactions. A hyd rostatic 
lower boundary condition is used (Zingale et al. 2002 1, imple- 
mented by Townsley and included in FLASH 33. ITie plane- 
parallel calculation is performed with a uniform 1 .22 x lO^cm 
grid of cells with large lateral width of 5 x lO^cm. We do not 
initialize any convective motion, rather we just use the thermal 
profile implied by a constant entropy and abundance convec- 
tive zone. 

We numerically initiate the detonation by raising the tem- 
perature in the convective zone to Tg = 3 in a circular region 
of radius 10^ cm, centered 1.2 x 10^ cm above the base of the 
convective layer This does not directly ignite a detonation, 
but rather creates a shock wave that traverses the cold layer; 
reflects from the inert core, and emerges out as a detonation 
wave|^ The resulting detonation front is shown in Figure [2j 

* http://flash.uchicago.edu 

^ This brings up a shortcoming of this calculation, which is that this inter- 
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which shows a snapshot of the detonation propagation 0.5 
seconds after the detonation was ignited at x = 4.5 x 10^ cm. 
The large lateral width of our simulation allows for the deto- 
nation wave to propagate far from the initiation site, demon- 
strating a steady-state form of much smaller spatial dimen- 
sions (< 10^ cm). The resulting speed of the detonation front 
is shown in Figure [3] by the points labelled "Blowout" and 
reaches 8.7x10** cm s"' for several tenths of a second as it 
crosses our large domain. 

The resulting composition profiles are shown in Figure |2] 
The burning is far from complete, with the "^He abundance far 
behind the detonation front only falling to « 0.5. Scarcely any 
""^Ni is produced, while more significant amounts of slightly 
lower Z elements are made, including '^'^Cr (bottom middle 
panel) and ''^Ti, but relatively little ^^Fe (not shown). The 
nucleosynthesis in the hot, carbon-enriched layer gives a sim- 
ilar amount of unburned helium, producing generally lower 
Z products such as '**^Ca (lower left panel), and appears to 
be a detonation wave partially supported by the overpressure 
from the underlying detonation. We performed a similar cal- 
culation in which the overlying layer had no carbon, and the 
detonation speed was slightly lower, 8.5 x 10^ cm s"\ but all 
other features were qualitatively similar The detonation front 
(vertically extending to 0.5 x 10** cm) is cellular, with the 
strongest shock-fronts propagating at an angle with respect to 
the overall leftward velocity. These cross-propagating shocks 
traverse both the cold pure helium and the carbon enriched 
layer above. 

We hypothesize that the surface blowout has induced a 
eigenvalue -type (sometimes called "pathological") detonation 
solution (Fickett & Davis 1979 1 in which the sonic point, 
where the flow becomes supersonic with respect to the shock 
front, is inside the reaction zone. This truncates the burning 
while still preserving a self-propagating detonation structure 
between the shock front and the sonic point (actually sonic lo- 
cus in multiple dimensions). This is surprising, as discussed 
below in Section[3] a detonation in uniform pure helium at the 
density at the base of the helium layer has a ZND length (to 
reach pure ^*Ni) much longer than our grid. 

3. LARGE AND TRANSIENT DETONATION STRUCTURES 

As we have shown (see also Sim et al. 2012), the low he- 
lium densities near the WD surface allow matter to blowout 
on a timescale comparable to the burning time. This is in 
contrast with earlier work done on Helium detonations in 
denser environments. For example, ,Khokhlov, (1989) per- 
formed ZND calculations for pure helium at ps = 50, find- 
ing complete burning by the time (and distance, 10^ cm) the 
flow behind the detonation front reaches the sonic point and a 
speed, 1.5 X lO** cm s"', consistent with the Chapman- Jouget 
(CJ) value for complete burning. However, at lower densities, 
the detonation structures become much larger than the verti- 
cal (and even ho rizontal) scale of the outer layers of a WD 
(Timmes & Niemeyer 2000), requiring a new approach. It 
also raises the possibility that a propagating detonation is not 
even possible in these thin helium layers, a quenching out- 
come well known in chemical detonations. 

We begin by investigating how marginally ignited detona- 
tions strengthen towards the steady-state CJ solution and how 

action led to immediate ignition of C/O when we allowed for carbon fusion. 
This is an artifact of our initiation scheme, and led us to use ^''Mg for the 
underlying layer in this initial study. 



this depends on composition. At low densities, these solu- 
tions will extend to physical dimensions much larger than can 
be accommodated on a WD, but will enable comparisons to 
ZND calculations that help us understand the likely effects of 
a finite geometry. In order to facilitate direct comparison with 
ZND calculations, we work in a 1 -dimensional plane parallel 
geometry in a medium of uniform density, temperature (typ- 
ically at = 0.1) and composition, sometimes including a 
small carbon fraction. The detonation is ignited by setting 
the temperature close to the boundary at = 3 with a dis- 
continuous fall to the ambient temperature. We usually use 
the smallest hot region that leads to a propagating detonation, 
decrementing a factor of 10 at a time. We perform these sim- 
ulations at two different resolutions, a coarse resolution to al- 
low the detonation to reach steady state, and a fine resolution 
which resolves all the burning stages. Adaptive mesh refine- 
ment is used in all cases, so that the resolution is actually the 
minimum cell size. 

We start by comparing to Khokhlov (1988 1989)'s calcu- 
lations at ps = 50 and Ty = 0.2 for the initial state. The speed 
at which the shock front propagates as the 10^ cm resolution 
simulation progresses is shown by the black curve in Figure 
[3] This detonation is of the CJ-type and the CJ speed (as de- 
termined by a calculation containing only the 13 a-chain el- 
ements) is shown by the horizontal dashed line at a speed of 
1.54 X 10^ cm s"'. Th is is also the spee d obtained by inter- 
polating in Table IV in [Khokhlov ( [1988 1, and we confirm his 
spatial profiles. Though it only takes 0.01s to reach 80% of 
the CJ speed, the full speed is not attained for about one sec- 
ond. 

We now consider a case with ps = 5, similar to the surface 
detonation calculation presented above. The evolution of the 
shock velocity (solid red) is shown in Figure [3] and is qual- 
itatively similar to the higher density result. Note that there 
are two red curves, one at a coarse (4.8 x 10^ cm) resolution 
and one much finer resolution (1000 cm). We performed the 
fine resolution simulations at a variety of ignition sizes, find- 
ing that ignition sizes of 10^ cm and below fail to produce 
a propagating detonation. We therefore use an ignition size 
of 10'' cm for all the fine simulations shown here. Similar to 
the higher density result, this calculation asymptotes to the CJ 
speed, but only after 1000 seconds and a propagation length 
w 10'^ cm, larger than the whole WD. 

We also expl ored the effects of a mod est fraction of '^C in 
the fuel, which [Shen & Bildsten ( 2009 1 suggested may pro- 
vide a faster reaction channel via '^C(a,7)'^0. We performed 
simulations with X\2 = 0.05 and 0.1, shown by the blue and 
magenta curves in Figure [3] Before about 10"^ s, the shock 
velocity is that of the leading shock from our ignition, not yet 
a detonation. The region between the contact discontinuity 
at the edge of our ignition region and this leading initiation 
shock is effectivly preconditioned by the shock propagation, 
and runs away to ignite the detonation. The spike in propaga- 
tion speed at a fewx 10"^ s is when this newly-formed deto- 
nation first passes the leading shock. It leaves this pre-heated 
region with an enhanced speed due to the preheating, falling 
quickly as the detonation propagates into now unheated fuel. 
Due to the presence of carbon, there is a continued enhance- 
ment of the propagation speed up to about 0. 1 s, potentially 
relevant if carbon is present in the accumulated cold He layer. 
As the detonation continues to strengthen, the speed becomes 
dominated by the consumption of "*He, and asymptotes toward 
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Fig. 2. — Steady-state temperature and abundance (in mass fractions) structure of a surface detonation propagating from right to left in the accreted helium 
layer of a WD. Prior to ignition, a hot, actively burning convection zone (0.25-1.4 on the vertical (radial) axis) enriched in '~C ashes sits above a cold helium 
layer (0.0-0.25) atop an inert core. Density contours are shown at 10^ lO", lO' and lO*" g cm . This snapshot is from 0.5 s after the detonation was ignited at 
X = 4.5 X 10^ cm. Only a portion of the computational domain near the detonation front is shown. 



the CJ solution. 

To investigate the cellular structure of these detonation 
fronts, we also performed a 2-dimensional calculation. The 
configuration was plane-parallel in a uniform medium, but 
now with random 1 % density perturbations near the ignition 
front to seed the cellular instability. The domain width trans- 
verse to the detonation front was 10^ cm with periodic bound- 
ary conditions, adequate to accommodate several detonation 
cell widths, though the cell width appeared to grow with time 
as the detonation strengthened. The resolution for the 2-d sim- 
ulation was about 5000 cm, enough to marginally resolve the 
fastest stage, helium capture on produced carbon. Even so, the 
evolution of the detonation speed (green curve in Figure [3]l 
agrees quite well with the 1-d simulations, giving us confi- 
dence that our 1-d simulations are robust enough to compare 
to the inherently multi-dimensional blowout case. 

This 2-dimensional uniform density simulation also allows 
us to study resolution effects in a more controlled context. 
By comparing a calculation at 7.8 x 10"* cm resolution to the 
one at 5000 cm, it appears that size of the cellular structure is 
influenced by resolution, with larger distance between cross- 
propagating shock fronts at coarser resolutions. However this 
does not have a significant impact on the resolved nucleosyn- 
thetic structures such as for the higher-Z elements (i.e. Ca, Ti, 
Cr, Fe and Ni). 

We find that weakly ignited (initially underdriven) detona- 



tions in helium strengthen toward the steady-state solution 
over time, but if the overall detonation structure - that is the 
thermodynamic and compositional profiles in space - is large, 
it will take a significant run of space for the detonation to take 
on this steady state structure. The incomplete burning struc- 
ture created by the blowout near the WD surface is similar to 
the early phases of these strengthening detonations, providing 
a valuable comparison. 

4. VERIFICATION BETWEEN ZND CALCULATIONS AND FLASH 
SIMULATIONS 

We have established that the helium detonation blowout at 
low densities reduces the available burning time, inhibiting 
the complete burning of helium and highlighting the impor- 
tance of accurately calculating the production of intermediate 
mass products. We begin by exploring how the details of the 
reaction network and calculation method impact the nucle- 
osynthesis. So as to save computati onal time, the FL ASH cal- 
culation uses the "approxl3" (Timm es et al.|20"00a| l reaction 
network, a simplified 13 nuclide network with steady-state ap- 
proximations of heavier p-rich nuclei along (a,p)(p,j) chan- 
nels and their inverses. We check whether this approximation 
leads to any substantial nucleosynthetic differences by com- 
paring to ID ZND calculations performed with a 200 nuclide 
nuclear network. 

In a standard ZND computation, the eigenvalue detonation 
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Fig. 3. — Detonation shock velocity as a function of time for detonations 
in lielium. Two densities for uniform fuel are shown, 5 X 10^ g cm"' (solid 
black) and 5 X lO' g cm"' (ID: solid red, 2D: green), and dashed lines show 
the corresponding steady state CJ speeds. The speed for material containing 
some '^C in fuel, Xjt = 0.05 (blue) and X[2 =0.1 (purple) is initially higher 
but asymptotes to similar behavior The propagation speed of the detona- 
tion wave moving through the helium surface layer on a WD (thick black) 
stabilizes well below the CJ speed. 

speed (or, equivalently, shock strength) is obtained when the 
position in the flow where the heat release function (see e.g. 
|Khokhlov|1 989 ) is zero coincides with the flow velocity being 
sonic with respect to the shock front. This gives a physically 
consistent transition from the detonation structure to the fol- 
lowing flow beyond the sonic locus. However, we cannot per- 
form a standard ZND computation for anything but the true 
steady state, which in this case is a far stronger detonation 
than either the early transient or blowout cases. Instead, we 
perform the first portion of a ZND integration, starting from 
the shock, for the known shock strength drawn from simula- 
tions and truncate our integration near the sonic point where 
it would become singular 

The first step is to examine the nucleosynthetic differences 
between a ZND calculation with approxlS and the 200 nu- 
clide nuclear network at the same conditions (ps = 5 and 
= 1). We choose to initialize the comparison ZND in- 
tegrations with the detonation velocity achieved at 0.2 sec- 
onds in the high resolution uniform medium 2D FLASH 
run, 1.09 x lO** cm s"'. This velocity is still lower than the 
steady-state CJ detonation velocity, therefore the ZND inte- 
gration will terminate when the downstream velocity in the 
shock frame is equal to the sound speed, here at a distance 
of 4 X 10^ cm. The resulting nuclide abundance distributions 
are shown in FigureH] showing that approxl3 provides a suit- 
able reproduction ofthe results obtained from the 200 nuclide 
network. '^"^Ti is overproduced by about a factor of two at in- 
termediate length scales, indicating that this nuclide is likely 
slightly overpredicted in our yields. 

A next check is to verify that the intermediate element pro- 
duction is consistent between the uniform density simulations 
and the ZND integration, both using the approxl3 network. 
We compare, in Figure [5] the simulation (dotted line for ID), 
which captures the begmning phase of the detonation while 
the shock is still strengthening, to the truncated ZND inte- 
gration (solid line for ZND). The excellent comparison shows 
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Fig. 4. — Comparison of post-shock nucleosynthesis of detonations in 
pure He computed in the ZND formalism using two different reaction net- 
works. The dotted lines are calculated using the aproxl3 reaction network 
and the solid lines from the 200-isotope network. The detonation veloc- 
ity was D = 1.09 X lO' cm/s for each run and the ambient conditions were 
po = 5 X 10^ g cm"' and To = 10^ K. 
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Fig. 5. — Comparison of laterally averaged 2-dimensional (dashed) and 1- 
dimensional (dotted) FLASH simulations of the He detonation structure after 
0.2 s of propagation. Also shown is a ZND calculation (solid), which uses 
the FLASH detonation speed at 0.2 s, 1 .09 X 10' cm s ' . Ambient conditions 
were po = 5 X 10^ g cm"' and Tq = 10^ K. Both FLASH and ZND are using 
the approxlB reaction network. 

that the truncated ZND is a suitable approach to understand- 
ing the strengthening detonations. Additionally, we need to 
relate the 2D simulations to those in ID. The laterally aver- 
aged 2D results are shown by the dashed line in Figure [5] 
remarkably similar to the ID results. Apparently the cellum' 
nature of the 2D detonation did not significantly change the 
burning length scale from the ID case. 
We have shown consistency in the detonation structure 
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computed by all of the following: (1) a truncated ZND cal- 
culation with a large network, (2) a truncated ZND calcula- 
tion with the aproxl3 network, (3) ID hydrodynamic simu- 
lations of the transient phase of a strengthening detonation 
that resolve all relevant reaction length scales, (4) 2D hydro- 
dynamic simulations that nearly resolve all relevant length 
scales. These consistencies give us confidence in the prop- 
agating solution revealed by the blowout simulations, but we 
have yet to use these truncated ZND calculations to fully di- 
agnose (and maybe predict) the blowout results (both nucle- 
osynthetic yields and speeds). 

5. CONCLUSIONS 

We have shown that a thin helium layer on the surface of a 
WD can host a steady-state, laterally propagating detonation 
at low densities (< 10^ g cm"^). Under these conditions, the 
detonation takes on the character of a eigenvalue-type deto- 
nation in which the sonic locus occurs well before complete 
helium burning is achieved. The stability and compact struc- 
ture of the resulting detonation is surprising because the size 
of the reaction zone of a steady-state detonation, which burns 
to completion, at these densities (3> 10^ cm) is larger than 
the layer thickness (^ 10^ cm) or even the star itself. Our 
2D reactive hydrodynamic simulations resolved the subsonic 
region behind the shock, exhibited the radial blowout of the 
hot material, and primarily yielded intermediate mass prod- 
ucts (i.e. ""^Ti and '^'^Cr) and relatively little ^^Ni. The in- 
complete burning had a spatial scale much smaller than the 
vertical thickness of the helium layer, and a detonation speed 
of 0.9 X 10^ cm s"', lower than the Chapman- Jouget speed of 
complete burning, 1.5 x 10^ cm s"'. 

This lower detonation speed reduces the strength of the 
inwardly moving shock that may ignite a d etonation in the 
carbon-oxygen core an d cause a Type la SNe ( [Finket al.|2007[ 
[20T0l [Simet al.|2012| l. The slower speed moves the conver- 
gence pomFTortEese weak shocks further towards the sur- 
face, potentially making it harder for the carbon to detonate. 
If a carbon detonation still occurs, the reduction in high-Z ele- 
ments that we found for the helium layer would allow for col- 
ors and spectral features consi stent with normal Type la SNe 
over most of their evolution (|Kromer et al.|[20T0| |Woosley | 

larger helium layer 



|&Kasenl20lT] [S 

masses. Modest 



im et al.'2012) for even large 



layer 

carbon enrichment in the helium layer was 
found here to produce significant amounts of '^''Ca, which, 
if the underlying core detonated, would appear in spectra at 
high velocities with severa l possible sources of asymmetry, a 
common feature of S Nla (Ka sen et al.||2003[ Mazzali et al.| 
[2005] |Tanakaetal.|2008 ; Foley et al.|2012| i. It is certainly the 
case that the extreniely low-mass helium shell considered here 



would only effect the earliest time colors ( Woosley & Kasen 



2011 Piro 2012 1, highlighting the importance of obtaining 



If the carbon does not detonate, then the only observable 
signatur es would be from the ejected ashes of the helium det- 
onation ( [Bildsten et al.|2007[[Shen et al.| 2010), a ".la" super- 
novae. For this M = Mq WD, the gravitational binding energy 
of the material at the surface is 3 x lO'^ergs g"', and rough 
integrals over the 2D simulation yielded burning products of 
X4 = 0.66, Xn ~ 0.054, X40 « 0.066, X44 = 0.079, X48 ~ 0.053, 
X52 = 0.026, X56 = 0.0067, and a smattering of other elements. 
This incomplete burning releases w 4.0 x 10'^ erg g"', so that 
the ashes are unbound and reach infinity with an average ve- 
locity of 4500 km s~', a rather low velocity compared to nor- 
mal supemovae. The small ejecta mass (O.O28M0) will lead 
to a short duration (likely less than 2 days near peak light) 
and the small amount of long-lived radioactivities imply a low 
luminosity (<C lO'*^ erg s"'). This rather low energy release 
brings up the interesting possibility that, for thin shells, a det- 
onation might not unbind the surface layers, though the CO 
interior might still shock-ignite. 

There remain several important uncertainties. The greatest 
is related to the initiation of the detonation in the convective 



early observations of SNe la (INugent et al.|201 1 1 IFoley et al. 

m% - — 



helium burning layer on the WD surface (Shen et al. 2010 
[Woosley & Kasen 2011). We assumed here that a detonation 
was ignited and just considered its propagation and stability. 
We intend to first explore the dependence of the detonation 
speed on both the He layer thickness and density. Coupled 
with our truncated ZND calculations (see section|4]l, we hope 
that these insights will yield a direct approach to predicting 
the minimum possible layer that will allow for a propagating 
detonation. This more complete understanding of the mech- 
anism by which the blowout of the surface layer leads to an 
eigenvalue-type stable detonation solution should allow for 
predictions of nucleosynthetic outcomes across a broad range 
of parameter space without expensive multidimensional sim- 
ulations. 
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